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Compass-Heisenberg model on the square lattice 
— spin order and elementary excitations 
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- Quantum phase transitions 

- Quantized spin models, including quantum spin frustration 

- Quantum error correction and other methods for protection against decoherence 

Abstract. - We explore the physics of the anisotropic compass model under the influence of perturbing 
Heisenberg interactions and present the phase diagram with multiple quantum phase transitions. The macro- 
scopic ground state degeneracy of the compass model is lifted in the thermodynamic limit already by infinites- 
imal Heisenberg coupling, which selects different ground states with Z2 symmetry depending on the sign and 
size of the coupling constants — then low energy excitations are spin waves, while the compass states re- 
flecting columnar order are separated from them by a macroscopic gap. Nevertheless, nanoscale structures 
relevant for quantum computation purposes may be tuned such that the compass states are the lowest energy 
excitations, thereby avoiding decoherence, if a size criterion derived by us is fulfilled. 
Published in: EPL 91, 40005 (2010). 



Strongly coiTelated electrons in transition metal oxides lead 
to rich quantum physics controlled by spin and orbital superex- 
change interactions that are complex and often intrinsically 
frustrated due to competing interactions |[Tl|2l . A high frustrat- 
edness is realized in the orbital compass model |[3}{7|, resulting 
in a large degeneracy of ground states. In contrast to SU(2) 
Heisenberg interactions which are not frustrated and isotropic 
in spin space, compass interactions are locally Ising-like but 
the spin component involved in the interactions depends on the 
bond orientation. Recent interest in the compass model was 
triggered by the observations that it has an interdisciplinary 
character and plays an important role in a variety of contexts. 
It could: (i) serve as an effective model for protected qubits re- 
alized by Josephson arrays ||6|, or (ii) describe polar molecules 
in optical lattices and systems of trapped ions fSj. First experi- 
mental successes in constructing special networks of Josephson 
junctions guided by the compass model have already been re- 
ported m. 

Materials with large spin-orbit coupling may give rise to 
compass spin interactions for some lattice structures |TOli lead- 
ing either to the compass or the Kitaev honeycomb model 
ifTTl . Numerical studies 112111311 and the mean-field approach 
lfT4l suggest that when anisotropic interactions are varied 
through the isotropic point of the two-dimensional (2D) com- 
pass model, a quantum phase transition (QPT) between two 
different types of directional order occurs. Recently the ex- 



istence of this transition, similar to the one found in the ex- 
act solution of the one-dimensional (ID) compass model ITSl . 
was confirmed using projected entangled-pair state algorithm 
lfT6l . This implies that the symmetry is spontaneously broken 
at the compass point, and the spin orientation follows one of 
two equivalent interactions, as concluded recently within the 
multiscale entanglement renormalization ansatz (MERA) ifTTl . 

In this Letter we introduce a generalized 2D Compass- 
Heisenberg (CH) model and investigate to what extent the de- 
generate ground states of the quantum compass model are ro- 
bust with respect to perturbing interactions which may occur 
due to an imperfect design of the system. General perturbations 
could therefore prohibit the conservation of non-local quanti- 
ties characteristic of the compass model ll6l [T2lll9| . We assume 
that such interactions are of Heisenberg type, as suggested by 
possible solid state applications 110111811 . We have found that 
the compass ground state is fragile and an infinitesimal Heisen- 
berg coupling is sufficient to lift its semi-macroscopic (expo- 
nential in linear size) degeneracy and to generate magnetic 
long range order, either ferromagnetic (FM) or antiferromag- 
netic (AF) one, which privileges a pair of columnar compass 
states as the ground state, while the other compass states sur- 
vive as finite energy excitations. Indeed, the phase diagram of 
the CH model, shown in Fig.[T] appears to be very rich and ex- 
hibits different QPTs between various phases of Z2 symmetry 
II20I triggered via softening either of spin waves or of a semi- 
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Fig. 1: Phase diagram of tlie CH model Eq. (1) with AF compass 
coupling Jz > 1. Square and diamond mark the isotropic compass 
points Jx ~ ± Jz, where in each case four ordered phases meet. For a 
phase $Q = Gz,C'2,Fx,- ■ ■,^ denotes the spin order depicted in the 
corresponding inset (i.e., F for a FM phase, G or C for AF phases) 
and a = x,y,z the spin orientation. The QPT between Fx and C'^ 
phases (solid line) is affected by quantum corrections to the classical 
transition (dashed line), see text. 



macroscopic number of quantum states. 

We consider a CH model of spins 1/2 on the square lattice 

i,3 

with {crfj] being Pauli matrices, Sij = {(7^?^, cr,?^^-, crf^}, and 
two types of nearest neighbour interactions along bonds in the 
2D lattice; {i) the frustrated compass interactions {Jx, Jr} 
which couple {erf } components along a-oriented bonds (the 
two axes are labelled a = x, z), and (m) Heisenberg interac- 
tions of amplitude /. We will use hereafter (p and Jc to define 
interaction parameters such that 



J, = Jc COS ( 



Jx = Jc sin ( 



(2) 



with Jc = 1 serving as a unit. Numerical results are obtained 
with Lanczos exact diagonalizations for periodic 2D rectan- 
gular clusters with even number of sites of up to A^ = 36 
(TV = LxXLz with L^ and Lz longitudinal dimensions, except 
for A^ = 18 and TV = 32 clusters which are rotated, see Ref. 
11211 ). Lattice symmetries (translations and C2v point group) 
and the <t^ — > — tr^ symmetry (for any r = {i,j)) are used to 
reduce the size of the Hilbert space; note that in contrast to the 
Heisenberg model ct^qj = J2i'^r^^ not conserved, but only the 
parity P = (-l)'^t°t/2. 

We first examine the effect of Heisenberg interactions J > 
favouring G, -type (Neel) AF order (Fig. [T]l. For AF anisotropic 
compass couplings J^ > \Jx\ spins orient along the z axis and 
form AF columnar states that are subsequently linked by posi- 
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Fig. 2: Spin structure factors (a) S^, and (b) 5^, for the A'^ — 36-site 
periodic cluster (L^ = Lz = 6) with isotropic AF compass couplings 
(Jx ~ Jz) for selected k — {kx, kz), see inset in (b), and different 
I/Jc- (c): Order parameter AIgz for cluster sizes TV = 16, 18,32 
for increasing I / Jcoi(L) ([5j. Solid (dotted) lines for an anisotropy 
parameter (jf>i — tajC^{Jx/ Jz) = 37r/20 and the isotropic compass 
model (ISO) (t> = 7r/4. 



tive I into the G^-AF structure. The spin structure factor 
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(normalized so that S^ = Jk,Q with Q = (tt, tt) for a Neel 
state) probes the onset of 2D spin order expected for increasing 
I/Jc, see Fig. 12 a). The compass regime with columns uncor- 
rected between one another, indicated by S^ being maximal at 
kz = TT and independent of kx 1221 . changes surprisingly fast 
into the Gz-AF order — the expected peak at k = Q grows 
rapidly with increasing I/Jc and is already of the order of unity 
for I/Jc — 0.01 ! In contrast, the peaks at k = Q in spin 
structure factors S^ [see Fig. I2b)] and S^ (not shown) grow 
slowly with I/Jc and are much smaller than Sq in the range 
of < I/Jc < 1, which reflects the spin anisotropy and Z2 
symmetry of the ordered Gz phase; eventually this anisotropy 
weakens for increasing I/Jc- 

The onset of the Gz-AF phase for I ^ Jc can be un- 
derstood using perturbation theory for small both Jx/ Jz and 
J/J2. The unperturbed Hamiltonian (ID Jz couplings) selects 
a manifold of 2^" column-ordered ground states, each column 
i = 1,- ■ ■ ,Lx possessing a 5 = 1/2 degree of freedom is rep- 
resented here by r^ = ±1, depending on the orientation of a 
reference spin in the column. 

In L'th order two neighbouring columns are flipped via the 
horizontal compass couplings 



3 

with a coupling constant 

Jcoi(L) = 8JzCl(J,/8Jz)^, 



(4) 



(5) 
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which gets exponentially small with increasing column length 
L = Lz for finite systems considered here. Here cl is a con- 
stant increasing with L (cl — 20, 252, 3432 for L = 4, 
6, 8 lfT2l ). Jcoi{L) vanishes in the thermodynamic limit for 
Jx/Jz < 1. explaining the (at least 2^^-fold) ground state de- 
generacy of the compass model. In contrast, Heisenberg in- 
teractions couple the {cf ,, crf^^j^} components in neighbour- 
ing columns, and this results in a first order perturbative cou- 
pling LzlTjTj_^_i which favours an AF arrangement between 
the columns j and (j + 1). This second term obviously dom- 
inates the first one ( Jcoi) in the thermodynamic limit as soon 
as / ^ 0, which thus ensures the onset of the Gz-AF long 
range order This is indeed confirmed by the evolution of the 
order parameter Mgz = Sf^ ^■. as a function of I / Jco\ shown 
in Fig. Oc) — while Jcoi varies over several orders of magni- 
tude between clusters considered, the values of Mgz are almost 
identical for a given ratio Jx/Jz- 

The isotropic ( J^, = J^) case is specific, as Sq^ saturates then 
close to ^ (rather than close to 1) for Jcoi <^ I < Jc, and so 
does Sq. In this case the compass interactions favour a mani- 
fold of (2^^ + 2^= )-degenerate ground states (in the thermody- 
namic limit), among which row-ordered states with spins along 
X (along with column-ordered states considered above); the 
perturbing Heisenberg interactions select in this manifold the 
four Neel states with ordered spin components either a^ (Gz 
phase favoured for Jz > Jx) or a^ {Gx phase for J^ > Jz)- 
This symmetry breaking was recently studied in the compass 
model (at 7 = 0) using the MERA lITTl . 

The perturbative treatment allows us to identify various or- 
dered phases which develop from the 7 = line in the phase 
diagram of Fig. [T] For instance, the C'z phase (with order pa- 
rameter Mc = SL ^0 is found when the dominant compass 
coupling Jz > favours AF-ordered columns while a small 
7 < couples neighbouring columns ferromagnetically. Thus, 
the compass model (7 = line) is at a transition between two 
phases, and four different phases meet at the isotropic points 
Jx = iJz (square and diamond in Fig. [T]), showing that both 
QPTs in the compass model evolve into transition lines between 
phases with different types of spin order 

Investigating the phase diagram of the CH model (Fig. [T]), 
perturbative approaches cannot be applied when Heisenberg 
interactions are comparable to dominant compass terms, espe- 
cially when they frustrate each other, e.g. for | J^l < J^, 7 < 
and 17/^^1 = 0(1). For large negative Heisenberg couplings 
FM phases are favoured — with spins being perpendicular to 
the z axis to avoid AF compass couplings. Depending on the 
sign of Jx, the most favourable spin orientation is either along 
X (Fx phase with order parameter Mpx = 5'fg q-. for Jx < 0) 
or along y {Fy phase with Mpy = SL ^x for Jx > 0), not 
found within the compass ground states manifold. By compar- 
ing the classical energies (per site) of different phases, we de- 
termined the classical transition lines in Fig. [T] For instance, 
the phase transitions between the C^ phase and both above 
FM phases (7^^; and 7^^,), were found using: Eq{G'z) = —Jz, 
Eo{Fx) ^ Jx + 27, and Eo{Fy) = 27. 

Strikingly, the phase diagram of the quantum model (Fig. [U 
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Fig. 3: Order parameters M^ of the C'- phase (AJc), and (a) of the 
Fy phase with 4> = vr/lO; (b) of the F^ phase with (j) = 377r/20. 
Ic indicates the transition between C'i and either (a) Fy phase, or (b) 
Fx phase — in (b) its classical counterpart 7j' is also shown, (c): 
Size-dependence of Ale for values I /Jx close to the transition dis- 
played in (a), (d): Momentum-resolved lowest energy levels per site 
(-E„(k) is the (n -I- l)*** state of momentum k) on a A'^ = 32-site 
cluster, across the transition displayed in (b). The ground state energy 
£■0(0, 0) shows avoiding crossing. 



differs very little from that found by comparing classical ener- 
gies. First, some transition lines are not modified by quantum 
fluctuations. This occurs due to additional symmetries at the 
classical transition line. An example is the case of J^: = J^ > 
and 7 > discussed previously, where the Hamiltonian is in- 
variant under a 7r/2 rotation of both spins and lattice along the 
y axis. The first order transition point between two symmetry- 
broken compass phases, found at Jz = Jx in the compass 
model lfT6l (square in Fig. [Hi, extends to a more conven- 
tional (still first order) transition line between the fully (clas- 
sically) ordered Gz and Gx phases, stable for either Jz > Jx 
or Jz < Jx- 

Similarly, on the line defined by 7 = —Jz/2 and Q < Jx < 
Jz (classical transition line between G'z and Fy phases), the 
Hamiltonian Eq. (1) has an additional U(l) symmetry ll23l . 
and therefore the corresponding QPT occurs necessarily along 
this line. We have checked by considering bond correlations 
(ctqCt^) for a ^ P that no other intermediate (e.g. spiral) 
ordered phase develops in this range of parameters. For both 
phases, the dependences of order parameters on I/Jz and on 
cluster size [shown for fixed (p = tt/IG in Fig. [3ja), with size- 
scalings of Mc in Fig.[3jc)] indicate a first order transition by 
the slope of order parameter which increases with A^ at the tran- 
sition, see Fig.[3ja). The scaling with increasing size indicates 
a discontinuity across the 7 = ~Jz/2 point, whereas both or- 
der parameters are equal at this point and extrapolate to a finite 
value in the thermodynamic limit. 

We emphasize that, contrary to what happens usually in first 
order transitions, no level crossing occurs in the ground state 
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for finite clusters (see also Refs. ifTSl and 1241 ). where the pre- 
ferred spin orientation evolves continuously from the z axis (in 
the C^ phase) to the y axis {Fy), but a level crossing occurs be- 
tween the two lowest excitations, as shown for the C, —F.^ tran- 
sition in Fig.|3jd). This is related to the distinct broken symme- 
tries in both phases in the thermodynamic limit: (i) translation 
symmetry breaking along z in the C^ phase results in a 2-fold 
degenerate ground state with k = (0, 0) and (0, tt) momenta, 
while {ii) the spontaneous breaking of Z2 symmetry in the Fy 
(or Fx) phase results in a 2-fold ground state degeneracy, both 
states with k = (0, 0) momentum [but with different parity 
P = (— l)'^t°t/^]. The crossing between the second (0,0) and 
the lowest (0,7r) states occurs exactly at / = —Jz/'^ in the 
C^ — Fy transition, for the above symmetry reasons. 

In contrast, some QPTs in the phase diagram of Fig. [1] e.g. 
the C^ — F^ transition (dashed line), do not present additional 
symmetries at the classical level and are affected by quantum 
fluctuations. The energies (per site) of C^ and F^ phase, eval- 
uated in second order perturbation theory, are: 



E{F,) = 



EoiC'z) - 
Eo{F,,) + 



Jt 



F' 



8 J, + 4/ 

'' z 

8 J^ + 12/ 



Jz~l 



(6) 



(7) 



They are equal on the solid line separating both phases in Fig.[T] 
Note that numerical estimations of the transition, with help of: 
(%) the crossing between both lowest excitations [Fig. Od)], and 
(ii) the increase/decrease of Mpx and Ale order parameters 
[Fig. [3b)] are not only consistent between themselves, but also 
in very good agreement with Eqs. (|6]l and (|7]i. 

Quantum fluctuations also shift somewhat the transition line 
between Gz and C^ phases (for J^ < —Jz)- Contrary to intu- 
ition, their contribution is here larger when neighbouring spins 
are aligned. Thus the FM phase is stabilized by them over the 
C'^ phase (at the F^-C'^ transition), and the Cx phase over the 
Gz phase (at the Gx-Gz transition). 

Next we analyze low-energy excitations for finite clusters 
which depend on the interaction parameters in a remarkable 
way. There are two fundamentally different types of excita- 
tions for the CH model Eq. ([T]i: {i) spin waves, i.e., coherent 
propagation of single spin flips, and {ii) column flips, where 
all spins of a column are flipped. First, we employ linear spin- 
wave (LSW) theory to estimate the dispersions of spin waves 
in various ordered phases with an adapted vacuum state (for e'- 
er G-like phases, a canonical transformation 1231 allows one to 
use a single type of bosons as for FM phases). With a Bogoli- 
ubov transformation, one finds spin- wave dispersion for the Gz 
phase (for \Jx\ < Jz and / > 0) 



Wk = 2^{2Jz +Jk+ 4/)2 - ( Jk + /k)2 



(8) 



with Jk = Jx cos kx and /k ~ 2/(cos kx+cos kz). The ground 
state is 2-fold degenerate in the thermodynamic limit with mo- 
menta k = (0, 0) and (tt, tt) as indicated in Fig. HI and there 
are two spin-wave branches with dispersion Wk and oJk+Q, re- 
spectively. The agreement between these two branches and the 
lowest spin-wave excitation energies obtained for finite clusters 
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Fig. 4: The two types of excitations in the Compass-Heisenberg 
model: spin-wave excitations (SWE) with the lowest energies for clus- 
ters of A^ = 16, 24, 32 sites (squares, circles and diamonds); the LSW 
dispersion Eq. ([8} shown by solid lines, with minima Sa at k = (0, 0) 
and k = (tt, tt)), see Eq. ^; and column-flip excitations (CFE) with 
energies Ec (triangles) for two values of Lz = 4, 6 (with L^ = 4), see 
Eq. dlOl ). Parameters I = O.lJc and (p — 37r/20 are corresponding to 
Gz phase in Fig. 1. 



with A^ < 32 sites (identified as excitations to the lowest states 
with P = —1) is satisfactory, see Fig. H) The minimum of 
spin-wave dispersion Ea or anisotropy energy, e.g. 



Fa = 2^{2Jz -Jx+ 4/)2 - {Jx + 4/)2 



(9) 



in the Gz phase with excitations displayed in Fig. |4] is gener- 
ally finite in an ordered phase, and vanishes along a transition 
line to another phase (except for the / = line). An exam- 
ple is the Gz — Gx transition line, where numerics and the 
LSW theory indicate a gapless spectrum, in contrast with the 
isotropic three-dimensional CH model which has a finite gap 
ifTSl . The transition lines characterized by spin-wave softening 
differ from the / = line — there, e.g. at the Gz — G'^ tran- 
sition, spin waves remain gapped and a softening of (2^ — 2) 
columnar excitations occurs. 

For finite clusters and for small enough I/Jz excitations to 
other compass states have lower energies than the spin waves. 
The simplest ones of these excitations, called column flip, con- 
sist of reversing all spins of a single column — these excitations 
are fundamentally different from spin waves which occur for 
Heisenberg systems with long range order. The resulting exci- 
tation energy is mainly due to Heisenberg couplings between 
each spin of the column and its two horizontal neighbours, and 
is given by 

E,^4Lz\Il (10) 

thus increasing linearly with column size. The smallness of the 
coupling Jcoi Eq. ^ between columns gives rather weak k- 
dependence of these excitations, as shown in Fig. |4]for L^ = 4 
and fixed Lx ~ 4. This coupling decreases with increasing 
column length, being Jcoi = 2.3 x 10^'^ for L^ = 4 and Jcoi = 
1.2 X 10"** for Lz = 6 (both with = 37r/20) — thus almost 



p-4 



Compass-Heisenberg model — spin order and elementary excitations 



no dispersion is seen for L^ = 6. These excitations are shifted 
above the spin-wave ones with increasing size L^ and play no 
role in the thermodynamic limit. 

A comparison of the column-flip excitations with spin waves 
is of importance for finite clusters since such column flips could 
be used for fault-tolerant quantum computing, where a qubit 
would be encoded in the orientation of a given spin column. A 
criterion for the use of such a device is that spin-wave excitation 
energies remain above those of column flips, i.e., Ea > Ec, 
which defines a column-flip regime. In Fig. |5] we present 
regimes of these distinct excitations obtained for three clusters 
in the case of AF interactions, with fixed number of L^: =4 
columns and increasing column length L, = 2,4,6. On the 
one hand, the size of columns and the anisotropy ratio Jx/ Jz 
determine the range of the column-flip regime, which requires 
large coupling anisotropy (small Jx/ Jz) and sufficiently short 
columns (since Ec oc L^). In the perturbative regime of small 
I / Jz and Jx/ Jz one finds the transition between these two dis- 
tinct regimes at 



LzI^^Jj^{Jz-Jx)+fL 



(11) 



The factor f^, accounting for finite-size corrections to the spin- 
wave dispersion, vanishes in the thermodynamic limit; but re- 
markably, even in the vicinity of the isotropic point {Jx < Jz), 
these corrections allow compass excitations to be the lowest 
ones for small enough Lzl-^l- On the other hand, Lz must be 
large enough for column-flip excitations of different columns 
to be sufficiently far from one another in Hilbert space, so 
that qubits remain well protected against local fluctuations and 
noise ||6|. Thus, for given values of {Jx,Jz} couplings the 
system size must correspond to a compromise between all con- 
straints above, in order to define correctly qubits with help of 
these column-flip excitations. 

In contrast to frequently proposed quantum computing 
schemes 161 1261 . here a column-flip switches between quasi- 
degenerate eigenstates (split by an energy oc Jcoi); thanks to 
their columnar character the fault tolerance of the compass 
model persists, unaltered by perturbating Heisenberg interac- 
tions — imperfections in switches should not harm information 
storage more than usual decoherence sources. 

One may also wonder whether the features discussed above 
are a consequence of the particular nature of perturbing interac- 
tions. We can for instance consider, instead of Heisenberg cou- 
plings, XY-type couplings — or rather XZ-type within our no- 
tation, i.e., introducing /i,^ (erf (T^+ erf (Tp for each pair of near- 
est neighbours, such that the a\ spin components do not appear 
in the Hamiltonian anymore. One finds that the order induced 
by perturbations at even infinitesimal Ixz persists since e.g. the 
neighbouring columns in column-ordered compass states are 
again coupled at first order in perturbation. The global phase 
diagram would mainly differ from that found for the CH model 
(presented in Fig. [U by the absence of Fy phase, and the tran- 
sition lines would not be affected by quantum fluctuations at 
all. Nevertheless, the column flips will also be the lowest en- 
ergy excitations for sufficiently small perturbation amplitude 
Ixz and system size, allowing here again for a possible design 
of a quantum computation scheme. 
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Fig. 5: Two distinct regimes of low energy excitations for J^ < Jz 
1251 : column-flip regime with Ea > Ec (lower part) and spin- wave 
regime with Ea < Ec (upper part), as obtained for = 37r/20 and 
for Lx X L, clusters with: Lj^ = 4 and Lz = 2, 4, 6. For each system 
size, Ea and Ec are calculated separately. 



Summarizing, we have shown that the macroscopic 2^ 
ground state degeneracy of the anisotropic L x L compass 
model is lifted by infinitesimal Heisenberg interaction /. The 
Compass-Heisenberg model has a rich phase diagram — for 
small |/| long-range order develops from a pair of compass 
states, and the remaining compass states are split off by an en- 
ergy Ec ~ 4L|/|. Thus the spin waves are the lowest energy 
excitations in a large system. For nanoscale structures of length 
L, however, the sequence of excited states can be reversed, with 
quasi-degenerate compass states being pushed below the spin- 
wave excitations, provided the product L\I\ is small enough. 
In this way decoherence of column-flip excitations by decay 
via spin waves could be avoided in quantum computation ap- 
phcations. 
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